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Abstract 



In classical random matrix theory the Gaussian and chiral Gaussian random matrix 
F</ . . models with a source are realized as shifted mean Gaussian, and chiral Gaussian, 

random matrices with real (j3 — 1), complex (j3 = 2) and real quaternion (/3 = 4) 
elements. We use the Dyson Brownian motion model to give a meaning for general 
(3 > 0. In the Gaussian case a further construction valid for f3 > is given, as 
the eigenvalue PDF of a recursively defined random matrix ensemble. In the case of 
real or complex elements, a combinatorial argument is used to compute the averaged 
characteristic polynomial. The resulting functional forms are shown to be a special 
cases of duality formulas due to Desrosiers. New derivations of the general case of 
Desrosiers' dualities are given. A soft edge scaling limit of the averaged characteristic 
polynomial is identified, and an explicit evaluation in terms of so-called incomplete 
;7_j ■ Airy functions is obtained. 

in 

£2 ■ 1 Introduction 

O 

(N 

It has long been observed that random matrix theory links aspects of probability, combina- 
torics, integrable systems and asymptotics. As a concrete example, to be further developed 
in the present paper, let us consider the mean value of the characteristic polynomial of a 
matrix from the Gaussian orthogonal ensemble (GOE). With X an N x N matrix of stan- 
dard Gaussians N[0, 1], the latter consists of symmetric matrices G = ^(X + X T ). It is 
simple to see [20] that the problem of computing (det(AIjv — C))c7 g GOE can be reduced to 
the combinatorial task of counting the number of permutations on N letters consisting of 
j 2-cycles and N — 2j fixed points, weighted by (— 1) J A J ' and summed over j, and as a 
consequence 

(det(A/jv - G)) GeGO E = 2~ N H N (\), (1.1) 



where H n {x) denotes the classical Hermite polynomial. In fact (jl.ip holds for any G = 
2 (A + X T ) with the elements of X all independent and having mean zero and standard 
deviation unity. 

On the other hand, from the viewpoint of integrable systems, we know that the RHS 
of (jl.ip satisfies a second order constant coefficient differential equation. This differential 
equation has a turning point at A = y/2N, and this leads to the asymptotic formula (see 

e.g. m) 

lim J- e ~ x2 / 2 (det(\I N -G)) GeGO E , _ , /n ,, =Ai(y), (1.2) 



7V->oo /^C 1 ) 



2Af+2- 1 /2AT-i/6 J/ 



where Ai(y) denotes the Airy function and 

C« =7r 1 /4 2 -iV/2 + l/4 (Ar!) l/2 7V -l/12_ 

As V2iV is the leading boundary of support for the eigenvalues, this scaling of the char- 
acteristic polynomial is thus seen to provide a statistical probe of the largest eigenvalues. 
It is furthermore true that the RHS of (jl.lh permits an integral form, which implies the 
so-called duality formula [3], [7] 

(det(AIjv - G)) G eGOE = ((A + ix) N ) xeN[0tl/ ^ , (1-3) 

and this allows for an alternative derivation of ()1.2j) . 

Our aim in this paper is to seek analogous inter-relations for so called matrix models 
with a source. Both Gaussian and chiral Gaussian classes of such matrices will be consid- 
ered. In the Gaussian case, a simple example of such a matrix is G + Go, where G € GOE 
and Go is the diagonal matrix diag [(//) r , (0) N ~ r ], (a) p denoting a repeated p times. The 
matrix Go is referred to as a rank r spiked matrix. Similarly, a simple example of a matrix 
model with a source in the chiral Gaussian class is 



[X )pxn Upxp 



+ 



Z + Z , (1.4) 



Vnxn (^-Ojraxp 

I )px« Upxp 

where the rectangular matrix X nxp consists of standard Gaussians and the entries of 
(Ao) n xp are equal to fi on the diagonal for the first r rows, and zero elsewhere. 

We begin in Section 2 by using combinatorial arguments to derive the generalization 
of (II. ID for both the Gaussian and chiral Gaussian ensembles with a source, in the case 
the elements are real (labelled (3 = 1), and the case the elements are complex (labelled 
(3 = 2). In Section 3 we use the Dyson Brownian motion model [14j to give meaning to 
the Gaussian ensemble with a source for general parameter /3 > 0, and in Section 4 we do 
similar for the chiral Gaussian ensemble with a source. In Section 5 duality formulas from 
[8j are introduced, which generalize (II. 3D to the case of the Gaussian and chiral Gaussian 
ensembles with a source. The main purpose of this section is to give new derivations of 
the duality formulas. In Section 6 we link the evaluations of the averaged characteristic 
polynomials from Section 2 with (a special case of) the duality formulas from Section 5, 
and we proceed to compute soft edge scaled limits, giving rise to a generalization of (II. 2p . 



2 The averaged characteristic polynomial for Gaussian and 
chiral Gaussian ensembles with a source: real or complex 
entries 

2.1 The shifted mean GOE and GUE 

Let H be a symmetric matrix, and require that all entries on and above the diagonal 
be independent with mean zero, and require that the variance of the elements above the 
diagonal equal 1/2. Let Hq be a fixed diagonal matrix, Hq = diag (si, S2, ■ • ■ , s n ). We seek 
the mean value of the characteristic polynomial of the matrix H + Hq. 

Proposition 1. Let H and Hq be as specified above. We have 

N 

(det(XI N -{H + Hq))) = ( n(A - s 3 + **))«■ (2-1) 

Proof. Let H = [/tij]i,j=i,...,jv- We have 

N 

det(A/jv ~{H + Hq)) = Y, < p ) \[{\p{1) ~ s mi) ~ h i,P(i))> ( 2 - 2 ) 

PeS N i=i 

where e(P) denotes the parity of P and 

A - Si, i = j 



1 0, i ^ j. 

The fact the entries of H each have mean zero tells us that after averaging the only terms 
in (12. 2D to give a non-zero contribution consist entirely of fixed points and 2-cycles. Each 
fixed point P(l) = I contibutes A — si, while each 2-cycle contributes 1/2. Let there be 
N — 2j fixed points, and j 2-cycles. The (weighted) number of ways to form j 2-cycles 
from 2j numbers is 

while the (weighted) number of ways to choose the fixed points is 

N-2j 

l<h<h<---<lN-2j<N fe=l 

Further, such permutations have parity (— I) 3 and so 

[N/2] ( N-2j 

(det((XI N -(H + HQ)))=Y(-^ J ^ 3{ -^ E II( A -^)- (2.3) 

j=0 •*' l<h<l 2 <-<l N -2j<N fe=l 



Using the facts that 






we see that (J2.3J) can equivalently be written in the form (|2.ip , D 

In the special case that Sj = for each j, we see that (|2.1|) reduces to (|1.3|) . 

Similar reasoning applies to the calculation of the average of the characteristic polyno- 
mial for the shifted mean Gaussian unitary ensemble (GUE). In this setting H is required 
to be an Hermitian matrix with all entries on and above the diagonal independent with 
mean zero, and with the variance of the sum of the real and imaginary parts of the off 
diagonal equalling 1/2. The averaged characteristic polynomial is again given by (I2.ip . 

2.2 The shifted mean chGOE and chGUE 

Consider the matrix structure (jl,4p and suppose the entries of X are all independent with 
mean zero and variance unity. Let Xq have entries on the diagonal equal to Sj for rows 
j = 1, . . . ,p, and all other entries equal to zero. In the case that Sj = for each j and thus 
Zq = 0(„ +p ) x ( n+p ) in (|1.4p we know that [20] 

(det(AI n+p - Z)) = (-1)VA"-^(A 2 ), (2-4) 

where L^ n {x) denotes the classical Laguerre polynomial. We seek to generalize this formula 
to the case of nonzero Sj . 

Proposition 2. Let X and Xq be as specified above, and let Z and Zq then be specified by 
(7gp . We have 

(det(XI n+p -(Z + Z )))= { —^ — - / t n -Pe- t oF 1 (n-p + l;-X 2 t)T\(t + S f)dt. 

T(n-p+l) Jo ^ 

(2.5) 

Proof. In the analogue of (|2.2p for the matrix structure (jl.4[) it is again true that after 
averaging the only permutations giving a non-zero contribution are those consisting entirely 
of fixed points and 2-cycles. Further, there must be at least n — p fixed points, due to (jl.4[) 
having n — p zero eigenvalues. 

The 2-cycles (j'ij'2) giving a non-zero contribution have j\ (j'2) corresponding to a row 
(column) of X. Each time (J1J2) is a diagonal entry Xj ip+ j 1 — Sj 1 there is a choice in the 
analogue of (12. 2p of averaging s 2 , (contribution s?) or of averaging x 2 i j 2 (contribution 
unity). In particular, if there are j 2-cycles (and thus n + p — 2j fixed points) we must cat- 
alogue the 2-cycles according to the occurrences of there being zero, one, two etc. diagonal 



3 

that the weighted contribution to the averaged characteristic polynomial is 



entries each with the choice of averaging s^ for the corresponding j. Doing this we can see 



+( E ^)(":j)glj)o-»)'+-+ E »?,-<}p.6) 

i<Ji<?2<p v y v y l<h<-<Jp<p 

The above consists of p + 1 terms. Labelling them r = 0, 1, . . . ,p, the corresponding 

combinatorial factor 

/ n — r\ (p — r N 



3 -rj \J -rj 

results from counting the number of ways of choosing j — r rows from n — r, independently 
choosing j — r columns from p — r, and further independently ordering the j — r rows. 

The averaged characteristic polynomial comes from summing (J2.6I) from j = to p. For 
given j, we replace j i— )■ p — j. Making use of the explicit power series form of the Laguerre 
polynomial, 



w-sC-;)^- 



3=0 

we read off that 

(det(A/ n+p - (Z + Z ))> = (-irX n -p{p\L;-P(X 2 ) + ( £ «?) (p ~ IV-L^X 2 ) 

i=i 

+( E 44)(p-2)^n( A2 ) + ---+ E 4-<}- (2-7) 

l<h<h<p l<h<-<l p <p 

To proceed further, we substitute for the Laguerre polynomials their integral represen- 
tation 



m\L a m (x) = — / e-'t^oF^a + 1; -xt) dt (2. 

r(a + l) Jo 



(a special case of (|5.1ip below). This allows the RHS of f|2.T[) to be rewritten 

\n-p A 2 /-oo , . v 

+* p - 2 ( E «)+■■■+ E »?,-<}*• 

i</i<; 2 <p i<Zi<---<z p <p 

Performing the sum gives (J2.5J) . D 

Note that setting s/ = 0, / = 1, . . . ,p in (|2.7p reclaims f|2.4j) . 



If we replace X in (jl.4p by X> , then we could consider the case of X having complex 
entries, with real and imaginary parts having mean zero and the sum of the variance of 
the real and imaginary parts equalling unity. The above working remains valid and thus 
the average characteristic polynomial is again given by (|2.5p . 

Generally, for the matrix structure (|1.4p with n > p there are n—p zero eigenvalues, and 
the eigenvalues come in ± pairs. Squaring the nonzero eigenvalues gives the eigenvalues 
of (X + Xq)^(X + Xo). We can therefore read off from (J2.5J) the corresponding averaged 
characteristic polynomial. 

Corollary 3. Let the entries of X have mean zero and expected value for their absolute 
value squared equal to unity. Let X$ be nonzero down it diagonal only, with Xq Xq = 
diag (//i, . . . , n p ). We have 

(det(XI p -(X+X ^(X+X ))) = > } / r-^-*o J P 1 (n-p+l;-At)TT(t+^)di. 

T{n - p + 1) J ax 

3 Gaussian matrix models with a source 

3.1 The three classical cases 

Let G be a random Hermitian matrix with real (/3 = 1), complex (/3 = 2) or real quaternion 
(/? = 4) entries. With the joint distribution of the elements of G proportional to e~^' 2 ) Tr G , 
such random matrices are said to form the Gaussian orthogonal (/3 = 1), unitary (/? = 2) 
and symplectic (/3 = 4) ensembles respectively. Note that in the cases /3 = 1 and 2 this 
definition coincides with the definition G = ^{X + X ) for X consisting of standard real 
(/? = 1) or complex (/3 = 2) entries. By definition the corresponding Gaussian matrix 
models with a source [7] are then specified as random matrices H = G + M, where M is 
Hermitian and has entries from the same number field as G. Note that the set of matrices 
H may equally as well be referred to as the shifted mean Gaussian ensembles. 

For future purposes it is convenient to replace M by e~ T M (r > 0), and to modify the 
variance in the joint distribution of the elements of G so that the latter is proportional 
to e~ ( - ,3 / 2 ) TrG /( 1_e T ). We then have that the joint distribution of the elements of H is 
proportional to 

e -(/3/2)Tr (H-e-Af) 2 /(l-e- 2r ). ( 3-1 ) 

We now take up the problem of determining the eigenvalue PDF of H J181 §11.1 & §13.4.3]. 
A fundamental result in relation to Hermitian matrices with real, complex or quater- 
nion real elements is that the measure for the change of variables corresponding to the 
diagonalization H = UDW, where D denotes the diagonal matrix formed by the eigen- 
values A := {\j}j=i,...,N and U the matrix of eigenvectors, is proportional to (see e.g. [181 



Prop 1.3.4]) 

N 

II lAfc-A/JJ^j (UUU). (3.2) 

l<j<k<N j=l 

Here (WdU) denotes the normalized Haar measure on real orthogonal (/3 = 1), unitary 
(/3 = 2) and unitary symplectic matrices (j3 = 4). Using (|3.2p we can then deduce the well 
known result that the eigenvalue PDF, G T (X;jT) say, with /2 = {/■*./ }j=i,...,jv denoting the 
eigenvalues of M, for the shifted mean Gaussian ensembles is proportional to 

Y[ \x k - xtfe-m^-fr^Uti j e m^ { uKWL) {uUu) (3 3) 

l<j<k<N 

Here 

A = diag(Ai,...,Ajv), L = diag(/ii,...,/xjv) (3.4) 

and 

/3 = /3/(2(l-e- 2 -)), t = e^. 
Denoting 

' ^(c/AC/t^^t^ =: o j*/P$-jl) (3.5) 

the (un-normalized) PDF (j3.3|) reads 



Xjfe-^rf-P^tiojiWfafrn. (3.6) 



ll |Afe - /\, ' -' -' ■ 'n-' |, 



As it stands (|3.6p only has meaning for /3 = 1,2 and 4. One concern of contemporary 
random matrix theory is to give constructions/ interpretations of such PDFs for general 
/3 > (see e.g. [E5])- In fact two approaches are possible — one as the solution of a Fokker- 
Planck equation, which we take up in the next section, and the other as the eigenvalue PDF 
of an ensemble of matrices defined recursively, which is to be considered in the Appendix. 
We remark that in the case of the general variance Wishart ensemble, the corresponding 
interpretations of the analogue of (|3.6f) . proposed as a candidate for the corresponding 
/3-generalized distribution in |28j, have recently been given in |19j . 



3.2 Brownian evolution interpretation 



We see from (J3JJ) that as r — > oo the element PDF of H reduces to that for matrices G. 
The corresponding eigenvalue PDF is then well known (see e.g. [i~8l Prop. 1.10.4]) to be 
proportional to e~@ u with 

1 N 
i=i ^<j<k<N 

A key observation, due to Dyson [15] , is that e~^ u can be interpretated as the Boltzmann 
factor for a classical gas of N particles on the real line with potential (|3.7p . interacting at 
inverse temperature (3. Generally the Brownian evolution of such a classical gas is specified 
by the Fokker-Planck equation 

dpr r ^ d fdU I d \ 

" £ ^' £ = 2Z^bA- + ^A7> (3 - 8) 

For /? = 1, 2 and 4 this equation was shown by Dyson [15] to specify the eigenvalue PDF 
for Hermitian random matrices with real (/3 = 1), complex (/3 = 2) and real quaternion 
(/3 = 4) elements, and independent real and imaginary parts themselves Brownian paths 
confined by a harmonic potential. Later (see e.g. [24] ) this specification was shown to 
be equivalent to choosing the elements of H according to the distribution (|3.ip . Note 
that the solution of ()3.8f) corresponding to (|3.ip has the delta function initial condition 
Yli = i 5(Ai — fii). Consequently, for general j3 > (|3.6p is then, up to a proportionality 
depending on r, the Green function solution of the Fokker-Planck equation (|3.8p . 

Development of the theory of the eigenfunctions and eigenvalues of (|3.8[) allows Jq 
in (|3.6p to be made more explicit. This requires introducing the symmetric Jack polyno- 
mials P K (z; a) where z = (zi, . . . , Zn), k denotes a partition of length less than or equal 
to N (we write £(k) < N), and a is a parameter. These polynomials are homogeneous, 
P K (cz; a) = c' K 'P K (z; a) for c a scalar, and can be uniquely specified in terms of the polyno- 
mial eigenfunctions of a certain differential operator (see e.g. [181 §12.6] for more details). 
For a = 1 the Jack polynomials coincide with the well known Schur polynomials, and for 
a = 2 (a = 1/2) they are the zonal polynomials for the real orthogonal (unitary symplec- 
tic) group. Introducing too the quantity d' K as in [HI eq. (12.60)], it is a known result that 
Fq in (|3.6p can be expanded 



CK 



making (|3.6p explicit for all /3 > 



> (x,y)-^a ^p K((1) ^ ;a) , l^-yj 



4 Chiral Gaussian random matrix models with a source 

The chiral Gaussian ensembles consist of the first of the matrices in (|1.4|) , with X T replaced 
by X* , and the entries of the n x p (n > p) matrix X independent standard Gaussian real 
(/3 = 1), complex {(3 = 2) and real quaternion (/3 = 4) random variables. The shifted mean 
version of this ensemble is then given by the sum of the matrices (|1.4|) . 

With M := Xq, as in (I3.ip it is convenient to introduce a parameter r, r > 0, such 
that Y = X + M is modified to have joint distribution of its elements proportional to 



-( / 9/2)Tr(y-e- T M)t(y-e- T M)/(l-e- 2T ) 



(4.1) 



We seek from this the eigenvalue PDF of (Y'Y) 1 ' 2 . Let the eigenvalues of the latter be 
denoted {yj}j=l,...,p- The strategy (see [181 §11.2.2]) is to introduce the singular value 
decomposition Y = UKV where U is an n x p unitary matrix, V is a p x p unitary matrix 
(both U and V have elements of the same type as Y) and A is an n xp matrix with diagonal 
entries y\,...,y p and all other entries zero. The measure for this change of variables in 
proportional to 



p v 



]>f +1 II bl-y^Rdy.iuUiDivUvi 

where a = n — p +1 — 2//3. 

Multiplying this and (|4.ip . then integrating out over U and V shows that the PDF for 
{yj}j=i j ... tP , G^. h (y;/2) say, is proportional to [18i eq. (11.105)] 

p 

n/3a+l TT I 2 2 1/3 -0T P ,v 2 -Bt 2 y p , a 2 

x ^ 1 (2//3) (/3n/2;{y]} i=li ..., p ;/3t{^ 2 }, =1 ,..., p ) (4.2) 

where J3, t are as in (13. 3j) . {/i 2 } are the eigenvalues of (M^M) 1 / 2 and 

^ 1 (2//3) (/3n/2;{y 2 } J=1 ,..., p ;/3t{ M 2 } J=1 ,..., p ) 

"(l/tdl7) /'( F t dF ) e ^iv(VAt [ /t i+itt/ Ayt)_ (43) 



If we introduce Xj = y 2 , rrij = ^i 2 (jf = 1, . . . ,p) as the eigenvalues of X^X and M^M 
respectively, then as a PDF in x := {a;j}j=i,...,p (|4.2p reads 

v 

Yl x T' 2 kfc -^l^e-^ E ^^-^ 2 ^=i m ^ -F 1 (2//3) (/3a/2 + /3(p- l)/2 + l;f j; ^m). 
jf=i i<j<fc<p 

(4.4) 



Analogous to the discussion in §3.2|, the fact that (|4.ip satisfies a multi-dimensional 
heat equation can be used [T7] to characterize (|4.2p as a solution of the Fokker-Planck 
equation (|3.8fl . with N t-> p, {Xj} i->- {y^} and 



^4£^^X>^ 2 - E logi^-^i, 



j=l 3=1 l<3<fe<P 

defined on the half line y > 0. Moreover, by developing the theory of the corresponding 
eigenfunctions and eigenvalues [3], we can show that for general a > 



-(a), ^ -* s^a^P K (x;a)P K {y;a) 
V [#V K P K ((l)*;a) 



(cf. (j!T5]> ) where 



'' r(c+(j-l)/a + Kj ) 



k ' 11 r(c+(j-l)/a) 



Thus (031) is explicit for all j3 > 0. 



5 Duality formulas 



The duality formula (jl.3p for the average of the characteristic polynomial in the GOE can 
be generalized to a duality formula for the n-th moment of the characteristic polynomial in 
the GOE [3j, [7J. In fact such formulas hold not only for the GOE, but their generalization 
known as the Gaussian /3-ensemble. To define the latter, let MEr n( w ( x )) refer to the 
matrix ensemble with eigenvalue PDF proportional to 

N 
Y[w( Xl ) H \x k -Xjf. (5.1) 

1=1 l<j<k<N 

The Gaussian /3-ensemble is then MEp t N(e~ x ). 

In terms of this notation, we have as a generalization of (jl.3p the duality formula |X8|, 
eq. (13.162)] 

N n 

3=1 ^ 3=1 

Note how the role of N and n is interchanged on the two sides, and that on the RHS 
(3 i— > 4//3. Desrosiers [8] (see also [7] in the case (3 = 2) has generalized (15. 2|) to the case 
of the Gaussian matrix model with a source. These formulas, together with derivations 
distinct from those given in [8], will be given in the next sections. 

10 



5.1 Gaussian matrix model with a source 



Let us denote the eigenvalue PDF (|3.6p for the Gaussian matrix model with a source in 
the case /3 = t = 1, by ME / g i Ar(e _A ;/T). In terms of this notation, the duality formula of 
[8] reads 

n N N n 

In writing this from the form given in [8] we have made use of the fact that the unspecified 
normalization in (|3.6[) is independent of //. 

The derivation of (j5.3|) given in [§] relies on special properties of the so called Dunkl 
transform. Here we will verify (|5.3p by using the generalized Hermite polynomials P« (x; a) 
[3]. The latter can be specified by the requirements that their expansion in terms of Jack 
polynomials has the structure 

pW (x; a) = P K (£; a) + ^J o K ^P^ (x; a) 

M:|mI<|k| 

for some coefficients a K/i , and that they exhibit the orthogonality 

Pf)(f;a)Pf'(x») oc<5 K , p . 

^ /ME 2/QiJV (e— 2 ) 

We require a number of special properties of the generalized Hermite polynomials, in 
particular integration and expansion formulas. The integration formulas required are [3] 

(^fe»>ME 2W (e^ ;l1 = ^(*5°)- ( 5 - 5 ) 

Similarly two expansion formulas are required. The first is [181 e Q- (13.130) upon the 
limiting procedure of Exercises 13.3 Q.5] 

N n 

n n^* - yi) = e^- 1 ) 1 " 1 ^^ «)« H " i/2 f ) (^ v«) (s.e) 

fc=n=i « 

where /i := (iV) n — k', with «;' the conjugate of k, while the second is its Jack polynomial 
analogue (see e.g. [El eq. (12.212)]) 

N n 

\\\[{x k + y l ) = Y J P^;a)P Pi {yMa). (5.7) 

k=ll=l K 



11 



Consider (|5.6p rewritten to read 

n N 

j=lk=l K 

Substituting in the LHS of (|5.3p and making use too of (|5.5p shows that the LHS of (|5.3p 
can be written 

Y J ^K^pP^a)Pl H \s ] a). (5.8) 

K 

Next we consider ([5.7p rewritten to read 

N n 

j = l fc = l K 



Substituting in the RHS of (15. 3p and making use of (|5.4p shows that the former is also 
given by (|5.8p . thus establishing the validity of (|5.3p . 



5.2 Chiral Gaussian matrix model with a source 

With f3 = t = 1, /3a/2 i— >• a, we will denote the transformed PDF for the chiral Gaussian 
matrix model with a source (|4.4p as ME^p(/e _I ;m). In terms of this notation, the duality 
formula of [8] reads 

n p 2 \ 

nTTfso + -s fc )) 
■ LJ - /ME fi „(i»e- I ;m) 

j=lfc=l M ^'^ ' ; 

n p 2 \ 

Ff 17(3; + ^m fc )) . (5.9) 

j=l fc=l ' 



A verification of (15. 9p analogous to that just given for (15. 3D is possible, again distinct 
from the derivation given in [8]. Our working makes use of the generalized Laguerre 
polynomials P^ (x;a;a) [3], specified by the requirements that their expansion in terms 
of Jack polynomials has the structure 

P^ L \x] a; a) = P K (x; a) + ^ a^P^x; a) 

p:\n\<\K\ 

for some coefficients a Kfl , and that they exhibit the orthogonality 

(p^\x-a;a)P {L \x-a;a)) ^5 K>p . 

\ ' /ME,/„ „(x a e-*) 
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With a = 2/(3, we seek to expand both sides of (|5.9p as a series in {P^ (—s; 1/a; a+a — 1)}. 
For this we require the Jack polynomial expansion formula (|5.7p . and the Laguerre 
analogue of the expansion formula (|5.6p , 



N n 



II 11^ - Hi) = ^(-l) M PfH%a;a)<x- M/ % L HVw?i V«;a(o + 1) - 1), (5-10) 

fc=U=l re 

which is derived in the same way as is indicated for (|5.6p . We require too the Laguerre 
analogues of the integration formulas (|5.4p and (|5.5|) [3j , 

(P«(^a)) M E a/aiW («ae-;£) = (-l) 1 " 1 ^-*; a; a) (5-11) 

(Pi L \x;a;a)) ME2/ ^ N(x a e - x . g) = P K (z;a). (5.12) 

Rewriting ([5.101) to read 

n p „ Oil 

II II ( a i + « Xfc ) = E(- 1 )' Ml ( «) ^ P - L) & 2 //3; 0a/2)pW (-s; (3/2; a + 2/(3 - 1), 

j=lfc=l ^ re P 

then substituting in the LHS of (|5.9p and performing the integration over if using (J5.12D 
shows that the LHS of (15.91) can be rewritten 



£(-l)Wp K (2m//3; 2/[3)P^ (-s; (3/2; a + 2/(3 - 1). (5.13) 

K 

To express the RHS in an analogous series, we rewrite (|5.7p to read 



n p 



[J [J(xj + 2m fc //3) = ^P K (2m//3;2//3)P M (f;/3/2). 

j=l fc=l K 

After performing the integration over x using (|5.1ip we reproduce (|5. 13[) . 

6 The averaged characteristic polynomial for Gaussian and 
chiral Gaussian ensembles with a source: general f3 > 

In this section the duality formulas ([5.3|) and ([5.9p will be used to generalize ([2.1|) and ([2.5|) 
to the case of the Gaussian and chiral Gaussian ensembles with a source for general (3 > 0. 
The resulting formulas will then be used to deduce scaled limits generalizing (II. 2p . 
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6.1 The Gaussian ensemble with a source 

Consider (15.31) with 



n = 1, s\ = A, 2/a = /3, x^-ix. 
Simple manipulation shows 

AT JV 

where ME^^e"^ 2 / 2 ; £) refers to the PDF jMJ with /3 = /3/2, i = 1, A = y, // = s. For 
/3 = 1 (/3 = 2) this is precisely (|2.ip in the case that i? has real (complex) off-diagonal 
elements. 

The RHS of (|6.ip is independent of /3. But in the case j3 = 2 and x i— >■ x/2 we know 
that the average characteristic polynomial gives the multiple Hermite polynomials of type 
II [5]. Thus this identification persists for general f3 > 0. 

In relation to the asymptotics we restrict attention to the case Xj = (j = r+1, . . . , N). 
Minor manipulation of (|6.ip shows 



n<* -»>>„. ..,.-,^ = (-»"l«' rM ^ <-**». 



where 



MEj, N (e-?» 2 /^ V ' 2 N 



/200 J 

e ,V4 + ^iV- r -Q (y _ a . ) |y (62) 

•ioo =1 ^ 7ri 

The function T^ r+1 '(u; {a*;}) is a particular incomplete multiple Hermite polynomial [9]. 
Introducing the scaled variables 

\ = V2N + -j=^, x k = ^N]2-—jJ^ (k = l,...,r) (6.3) 

the leading large N asymptotic form of the incomplete multiple Hermite polynomial can 
be read off from ^ Proof of Prop. 18]. This involves the incomplete multiple Airy function 

( _ ir +i Ai (r + D (X){sfc}) := j e -WV3n (w + Sfc > 

= ri(-7^+^) Ai ^' ^ 



dx 

fc=i 

where *4 is a simple contour, starting at e~ nl ' 3 oo, and finishing at e^^oo, following the 
corresponding rays asymptotically. As a consequence, we can generalize (jl.2p to the case 
of the Gaussian ensemble with a source. 

14 



Proposition 4. Let A and x k be given in terms of X and s k according to i6.3\) . Define 



c (r+l) _ 2 -(N-l)/2 N (N+l)/2-(r+l)/3 e -N/2_ 



We have 



~\ 2 /2 N 

Km ^—-(H(\-y k )) =(-iy+iAi(r+V(X,{s k }). (6.5) 

Ls N k=l 3,a/v , ; 

Making use of Stirling's formula to replace (TV!) 1 ' 2 in (jl.2p . and noting from (J6.4J) that 
-ArW(X, {s fc }) = Ai(X) shows that the r = case of dESD agrees with (TOD . 

6.2 The chiral Gaussian ensemble with a source 

With n = 1, si = —A, m i— >• f3rh/2, a i— >• (/3/2)(a + 1) — 1 and the change of variables 
2xfc//3 I—?- a?fc on the LHS, the duality formula (|5,9|) reads 

n(A-x fc )) 
7 ME« ( I (W2)(«+l)-l e -/3x/2 ;W 2) 
fc=l 



L T-7T / x a e" x 0^1(0 + 1; -Ax) TT(x + m fc )da;. 

r(a + i)7 ^ 



(6.6) 



We recognize the RHS of (|6.6p as the multiple Laguerre polynomial of type II [ITJ eq. (77)]. 
Also, we remark that this evaluation of the averaged characteristic polynomial, for a = n—p 
and /3 = 1 or 2 is in precise agreement with the result of Corollary [3l 

A systematic study of the asymptotics of the multiple Laguerre polynomials is yet to 
be undertaken. We will initiate this task in the simplest nontrivial case of 

7712 = • • • = 77lp = 0. (6-7) 

Proposition 5. Let J 6', 7| ) hold, and introduce X and s\ so that 

X = ip + 2a + 2 + 2(2p) 1/3 X, m x = p - {2p) 2/3 s x . (6.8) 



Also define D [ p > = pl(2p)^ 3 2- a . We h 



ave 



e -V2 / P v 

lim — -( r[(A-x fc )) = siAi(X) - Ai'(X). (6.9) 



Proof. Recalling (|2.8p we see that the RHS of (|6.6p can be written in terms of Laguerre 
polynomials to read 

(-l)V^(A) - m^-iy-'ip - ly.L^X). (6.10) 
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But it is a consequence of an asymptotic formula in Szego's book J27J p. 201] that for fixed 
k and with A as in 



/ Ik 

e- x '\-if +k v; )+k {\) = 2- a - 1 /3 p -i/3( Ai(x) _ __ A i'(x) + o( P - 2 / 3 ; 

Multiplying (|6.10p by e ' 2 , substituting for m\ according to (|6.8p and using this asymp- 
totic formula with k = and k = — 1 gives (|6,9p . □ 



Note from (|6.4p that the scaled limit (|6.9|) is the same as the scaled limit (|6.5|) in the 
case r = 1 -- an example of universality. More generally, choosing rrik = p — (2p) 2 ' 3 Sk, 
k = 1, . . . , r, and m^ = 0, fc = r + l,...,pwe would expect that for an appropriate Dp , 
that upon multiplication by e~ ' 2 /Dp , the large p limit of the averaged characteristic 
polynomial in (|6.6p will equal the RHS of (|6,5p . 



6.3 Related work and future problems 

In the Introduction it was remarked that, in the classical setting, one link with integrable 
systems is provided by the fact the Hermite polynomial in (jl.ip satisfies a differential 
equation. It is also true that the incomplete Hermite polynomial (|6.2p satisfies a linear 
differential equation — now of degree r + 2 (see e.g. [25, eq. (23.8.7)]). A future challenge 
is to find a derivation of this differential equation, working directly from the definition of 
the random matrix average in (|6.ip . Such a derivation is quite straightforward in the case 
of the average in (jl.ip [1], upon making use of the explicit form of the joint eigenvalue 
PDF, and is in fact perhaps the simplest result in the broader concern of finding (partial) 
differential equations associated with Selberg correlation functions [26], [22j . 

Another way to characterize the Hermite polynomial in (jl.ip is via its three term 
recurrence. The tridiagonal construction of the Gaussian /3-ensemble due to Dumitriu and 
Edelman leads most naturally to a direct demonstration that the characteristic polynomial 
on the LHS of (II. ip satisfies the same three term recurrence [13] . It is similarly the case that 
the multiple classical polynomials corresponding to the RHSs of (16. ip and (16, 6h satisfy linear 
recurrences [2] , thus suggesting the problem of a direct derivation of the same recurrences 
for the characteristic polynomials on the LHSs. For the shifted mean Gaussian ensembles 
of Section 12. 1[ with s r +i, . . . ,sjsf = 0, the r x r block tridiagonal form found recently by 
Bloemendal and Virag [6], may provide some insight. 

The duality formulas (|5.3p and (|5.9p give formulas well suited to asymptotic analysis, 
not just for the averaged characteristic polynomial, but for the average of an arbitrary (but 
fixed) number of characteristic polynomials and thus by coalescence for integer moments 
of the characteristic polynomials. Using the duality formula (15. 2p the soft edge scaling 
of the latter for the Gaussian /3-ensemble was given in terms of a /3-generalization of the 
Kontsevich matrix Airy function [10] , and for (3 = 1 , 2 and 4 (general /3 > 0) this was 
generalized in [8] ([32]) to the case with the variables not coalesced. Of course one would like 
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to compute the generalization of the soft edge scaling limit (|6.5p from a single characteristic 
polynomial to a product of characteristic polynomials. In has been announced by Desrosiers 
and Liu in |12] that they are preparing a paper tackling this problem. 
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A A recursive construction of the Gaussian /3-ensemble with 
a source 

The Gaussian and chiral Gaussian matrix models with a source are specified by the PDFs 
proportional to (|3.6p and (|4.4p respectively. We have seen that both can be interpreted for 
general /3 > as Green function solutions of particular Fokker-Planck equations describing 
the Brownian evolution of the underlying log- gases. Here we will specify the recursive 
construction of a random matrix ensemble having the eigenvalue PDF (|3.6p . While no such 
construction is known in relation to (|4.4|) . the PDF specifying the general variance Laguerre 
/3-ensemble can similarly be generated as the eigenvalue PDF of recursively defined random 
matrices [19j . 

Consider random matrices specified by the distribution on elements (|3.ip . but with the 
r-dependence scaled out so that it reads 

e -(/3/2)Tr(H-Af) 2 _ ( A _-q 

Let us take (3 = 1 (real entries) for definiteness, and furthermore write Y = Y/v to indicate 
that the matrix Y has size N x N. Then by first noting that (lA.ip can equivalently be 
written 

e -(f3/2)Tr(H N -D N )^ ( A _ 2 ) 

where D^ = diag M^ = diag (/j,i, . . . , /ijv)> as observed recently in [4] we have that 



ePBFH N+1 = ePDF 



X U X 12 ••• £17V+1 

xi2 A^ ••• 



n \( N ) 

xin+1 ... A^ 



(A.3) 



Here ePDF refers to the eigenvalue PDF, {A.- }j=i 1 ... 1 n denotes the eigenvalues of Hjq and 

i n = N[ m i,l], Xlj =^MV2} (j = 2,...,iV + l). 

In fact in (IA.3P the eigenvalue PDF of Hn + \ is fully determined by the distribution 

\ Xlj \ 2 = r[/3/2;l], (A.4) 
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where r[s;<r] refers to the gamma distribution (PDF proportional to x s ~ 1 e~ s ' a , s > 0). 
This remains true for j3 = 2 and 4 (the RHS of (|A.3|) should then be scaled by l/\//7, but 
this detail can be ignored). 

But (|A.4|) makes sense for general (3 > 0, and allows us to use (|A,3p to give a recursive 
generation of the eigenvalue PDF /3-ensemble coinciding with the eigenvalue PDF of the 
shifted mean Gaussian ensembles for f3 = 1, 2 and 4. The use of (JA.4J) to /3-generalize clas- 
sical random matrix ensembles goes back to |13|). while the idea of a recursive construction 
can be found in |21| . According to (j3.6[) the eigenvalue PDF of H^ is proportional to 

J] |A fe -A/e-^^/ 2 -^^?/ 2o ^) (X;/2) (A . 5) 

l<j<k<N 

for /3 = 1,2 and 4. We would like to show that the recursive construction (|A.3|) implies 
that (|A.5j) remains true for general /3 > 0. Our strategy will be to show that (1A.3P implies 
a recurrence formula in A~ for the eigenvalue PDF, and that (|A.5|) satisfies the recurrence. 

Proposition 6. Consider the random matrix on the RHS of iA.3\) . Let the eigenvalues be 

denoted {X, }j=i,...,N+i- With {X- }j=i,...,jv given, the eigenvalue PDF is proportional 
to 

e-CLf+M™* ) 2 -Ef =1 (Af ) f) e m + i (EfJl 1 *f +1) -£f =1 a< jv) ) 

y (iV+l) _^(A f +l)s 1 TV 7V+1 

Til (\W \ {N) )P-- 

lll<j<k<N\ A j ~ A k I j=l k=l 

subject to the interlacing 






A (»+d > A («) > A (» + i, > A («i > . . . > A <«) > A («+i). (A . 7) 

Proof. It is a simple exercise to show that the secular equation for the matrix on the 
RHS of (O is 

N 

= A-xn-^ — %-—, q j = \x hj+l \ 2 . (A.8) 

i=i A ~~ X j 
The distribution of the zeros of this random rational function in the case of x\\ fixed is 
given by [21] (Prop. 8 with n = N, s± = ■ ■ ■ = sn = /?/2, A i — ^ A — xn, a, t-t aj — Xu). 
Explicitly we have that in this setting the PDF for the zeros is supported on the region 
(|A.7|) and is subject to the further constraint 

N+l N 

E^ +1) = E^+-n. (A.9) 

i=i i=i 
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Subject to these requirements, its functional form is proportional to 

e - (Ef=V (Af +1) ? -Ef =1 (Af ' ?) e x\j2 

u<j<fc<JV+i^ A j A fc ! TT TT i \ W \ (JV+i) 1/3/2-1 



lll<j<fc<JV+lv A j A 

lli<j<fc<jv+i^ A j A fe J i=i k=i 



yUU^j-^t 2 - 1 - ( a - 10 ) 



Multiplying this by the PDF of ecu, expanding the exponent and making use of ()A.9|) gives 
(jMD - D 

Corollary 7. The joint PDF of {X- }j=i,...,tv and {A^- }j=i,...,jv is proportional to 

e -(EiL + 1 1 Of +1) ) 2 /2 e ^ +1 (Ef. + 1 1 Af +1) -Ef =1 Af ') "Q (A (AT+1) _ A (iV+l) ) 

i<j°<fc<iv+i 

AT iV+1 

>< n (A^-Af^nni^-^^i^-Vo^^;^) (a.h) 

l<j<k<N j=l k=l 



subject to the interlacing (A. 7). 

Proof. We simply multiply together (JA.6P and the PDF MEb n(^~ x ). □ 

To proceed further requires knowledge of (a special case of) the so-called Dixon- 
Anderson density (see e.g. |18l eq. (4.11)]). This is the conditional PDF for {A^- } given 

{A ■ } (note that (JA.6P is the other way around) specified by 

r({N + l)/3/2) ni<j<fc<Jv( A i ~ x k ) TT TT i\ (JV) \( jV + 1 ) i/3/2-i (A 19 n 
(r(B/2)) N + l n (\( N+ V \( N + 1 )\ H H ' j k ' ' { ' 

\ \H/ )) LLl<j<k<N+l\ A j - A fc ) j = l k=l 

and again subject to the interlacing (|A.7|) . Referring to this as DAj\r + i(/3/2), it is a known 
result that [23] (see also J2TJ Appendix C]) 

e „ +1 ES 1 xr 1} (e -- +1 EJL Af » oJ? /« (# ) ; ^ ))DANim 

= ^ 2//3) (A( iV+1 );^ iV + 1 )). (A.13) 

On the other hand, it follows from (jA.lip that the N i— >• N + 1 version of ()A.5p is also 
the PDF of {A,- } for general (3 > 0, if we can show that 

Af N+l 

^We^CEf^Ar^-Ef^A^) JJ (Af } - Af 5 ) JJ II |Af } - Af +1 V /2 ^ 

l<j<k<N j=l fe=l 

x ^ ) (A( Ar+1 ); M ( Ar + 1 )) 

oc n ( A f +1) - A f +1) ) /9 - 1 o^ //9) (^ +1) ;M (iV+1) ). (A.14) 

i<i<fe<Ar+i 
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But (JA.13P is precisely (|A.14|) . complete with normalization. It thus follows that (|A.5|) is 
indeed the eigenvalue PDF of Hat+i- 
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